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ABSTRACT 

Context. Remnant radio galaxies represent the final dying phase of radio galaxy evolution in which the jets are no longer active. 
Remnants are rare in flux-limited samples, comprising at most a few percent. As a result of their rarity and because they are difficult 
to identify, this dying phase remains poorly understood and the luminosity evolution is largely unconstrained. 

Aims. Here we present the discovery and detailed analysis of a large (700 kpc) remnant radio galaxy with a low surface brightness 
that has been identified in LOFAR images at 150 MHz. 

Methods. By combining LOFAR data with new follow-up Westerbork observations and archival data at higher frequencies, we 
investigated the source morphology and spectral properties from 116 to 4850 MHz. By modelling the radio spectrum, we probed 
characteristic timescales of the radio activity. 

Results. The source has a relatively smooth, diffuse, amorphous appearance together with a very weak central compact core that is 
associated with the host galaxy located at z=0.051. From our ageing and morphological analysis it is clear that the nuclear engine 
is currently switched off or, at most, active at a very low power state. We find that the source has remained visible in the remnant 
phase for about 60 Myr, significantly longer than its active phase of 15 Myr, despite being located outside a cluster. The host galaxy 
is currently interacting with another galaxy located at a projected separation of 15 kpc and a radial velocity offset of ~ 300 km s^L 
This interaction may have played a role in the triggering and/or shut-down of the radio jets. 

Conclusions. The spectral shape of this remnant radio galaxy differs from most of the previously identified remnant sources, which 
show steep or curved spectra at low to intermediate frequencies. Our results demonstrate that remnant radio galaxies can show a wide 
range of evolutionary paths and spectral properties. In light of this finding and in preparation for new-generation deep low-frequency 
surveys, we discuss the selection criteria to be used to select representative samples of these sources. 

Key words, galaxies : active - radio continuum : galaxy - galaxies : individual: J1828-)-4914 


1. Introduction 

Radio-loud active galactic nuclei (AGN) can have a substantial 
energetic impact on both their host galaxy and the surround¬ 
ing in tergalactic medium ([Briiggen & Kaiser|2002| |Croton et al.| 
2006[ Morganti et al.|2013| l. The jets of relativistic particles pro¬ 
duced by the AGN can reach distances of up to a few Mpc from 
the host galaxy ( |Willis et al.||1974{ |Machalski et al.|2008[ ). The 
active phase of the radio AGN can last several tens of Myr, af¬ 
ter which the nuclear activity stops and the source starts to fade 
away ( [Parma et al.|1999| l. In this fading period, often termed the 
remnant or dying phase, the core and the jets disappear and only 
the lobes remain visible, radiating their energy away ( jSlee et aL] 
|2001| l. When an accurate knowledge of the radio spectral charac¬ 
teristics is available, remnant radio galaxies provide a timing of 
the on-off cycle of the radio activity, which constrains the evolu¬ 
tionary history of the radio source ( |Kardashev| 1 962| [Pacholcz^ 
|1970[ jJaffe & Perola|[T973j l. The fate of the AGN remnant ra- 


* brienza@astron.nl 


dio plasma and the physical processes driving its evolution are 
widely relevant to studies of radio galaxy feedback. Moreover, 
a better knowledge of these objects could also give new insights 
into the formation and evolution of cluster AGN relics and radio 
phoenices ( jSlee et aL 200 1[ EnBlin & Bruggen|2002[|van Weerenj 


et al.|200^ de Gasperin et al.|2015|l. 


Statistical studies based on number densities of radio sources 
show that radio-loud AGN are active for a short interval of time 
with respect to the entire life of the host galaxy. The duty cy¬ 
cle of the radio-loud AGN strongly depends on the power of the 
radio source ( [Best et ^|2005| jShabala et al.||200^ . Therefore, 
we would expect to observe many sources showing signatures 
of remnant emission from past radio-loud AGN activity. But de¬ 
tections of such objects have always been elusive and rare. jGio^ 
vannini et al. ( 1988| l identified only a few percent of the radio 
sources in the B2 and 3C samples as remnant radio galaxies. The 
reason for this rarity is most likely connected to the rapid lumi¬ 
nosity evolution of the fading plasma ( Kaiser & Cotter|[2002| l. 
Synchrotron cooling and inverse Compton scattering of cosmic 
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microwave background photons cause a steepening of the spec¬ 
trum at high frequency that renders the source undetectable. In 
addition, adiabatic expansion of the lobes contributes further en¬ 
ergy losses in an energy-independent way, decreasing the flux 
densities across the whole spectrum. Although still poorly un¬ 
derstood, the timescales of these processes depend on the phys¬ 
ical conditions within the lobes and the ambient medium as the 
jets switch off. For example, if the lobes remain strongly over¬ 
pressured at the end of the active phase, adiabatic losses can be 
significant, unlike the case where the lobes have already reached 
pressure equilibrium with the surrounding gas. Moreover, pow¬ 
erful sources with high magnetic field values are subject to much 
higher radiative losses. 


Systematic studies of remnant radio galaxies as a popula¬ 
tion have always been hampered by the small number of sources 
available and by non-uniform selection effects that limit the 
study of their characteristics and evolution. The known objects 
cover a variety of morphologies, spectral characteristics, and 
physical conditions. Dying radio galaxies typically show relaxed 
morphologies without compact components and curved steep 
spectra Pamrozy et al. 2004[|Giacintucci et al.|2007[[Parma et'af] 


|2007||Murgia et al.|201 l| l even at frequencies < 1 GHz. In some 

cases, the fossil lobes look like bubbles of emitting plasma with 
a very irregular morphology (e.g. Giovannini et al.|[l991| |Har-| 


jris et al.|1993| jde Gasperin et al.|2014| l. Since radio jets can be 

episodic in nature, remnant plasma can also be associated with 
a pair of new-born jets ([Clarke & Burns] 199 1| [Jones & Prestonj 
[2001 [ [Saikia & Jamrozy|2009| l. These sources are often termed 

restarted radio sources. The most common example of these are 
the double double radio galaxies (DDRG, [Schoenmakers et ak] 
|2000| [Nandi & Saikia|2012| l that posess a pair of diffuse, steep- 
spectrum lobes which is well detached from the renewed radio 
activity. Other evidence of this occurrence is seen in the diffuse, 
extended emission surrounding some GHz-peaked spectrum and 
compact steep spectrum sources ([Baum et al.|1990[ Stanghellini 


[et al.||2005j [Shulevski et al.[|20i2[ l. We note that most of the 
remnant sources have been found in massive clusters in which 
the dense intergalactic medium prevents strong expansion losses 
( [Murgia et al.[[20lT| ). However, a few examples in low-density 
environments have a lso been identified pamrozy et al.|2004[ |de| 
Gasperin et al.|2014[ Hurley-Walker et al.|2015 1. 


To date, a unified physical framework in which to inter¬ 
pret the variety of individual studies mentioned above has not 
been provided. To constrain models describing the final stages 
of the life of radio galaxies, we must take a statistical approach 
( Kaiser & Cotter[[2002[ ) and, therefore, must expand the search 
for remnant radio galaxies. Greater statistics of their character¬ 
istics will enable an investigation of the physical processes driv¬ 
ing the remnant plasma evolution as a function of source power 
and environment. With this aim in mind, we plan to exploit 
upcoming deep and wide field radio surveys to build up com¬ 
plete samples of such sources. These new facilities are designed 
to significantly enhance the detection probability of these rare 
low-surface brightness objects. In particular, remnant emission 
is expected to be brighter at low radio frequencies and, there¬ 
fore, instruments such as the Low-frequency Array (LOFAR, 
[van Haarlem et ar]2013[l an d the Murchison Widefield Array 
(MWA, Lonsdale et al.|2009] ) will play a key role in the system¬ 
atic search for them. However, it is difficult to identify remnant 
radio galaxies, and simple selection criteria based on integrated 
spectra alone can result in a biased representation of the remnant 
population. For this reason, it is important to address the issues 
of remnant selection and associated selection effects. 
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A complementary approach to remnant radio galaxy selec¬ 
tion is radio morphology. We have preliminarily explored the 
available LOFAR fields to search for large radio sources with a 
low surface ratio without compact features, which may be in¬ 
dicative of remnant radio galaxies. Here we present one 700 kpc 
remnant radio galaxy at coordinates J2000.0 RA 18h 28m 20.4s 
Dec -i-49d 14m 43s (we call this blobl) that was discovered as 
part of this search, blobl was first identified on purely morpho¬ 
logical arguments based on its unusually large angular extension 
with amorphous shape, low surface-brightness, and lack of com¬ 
pact features. We investigate the physics of blobl in comparison 
with already known active and remnant radio galaxies with the 
aim to reconstruct the evolution history of the radio source from 
its birth to its death. This is also the approach that we intend 
to use for future sources. Furthermore, we discuss the selection 
criteria to be used for extracting complete samples of this pop¬ 
ulation of sources out of the upcoming deep low-frequency sur¬ 
veys. Such studies are an important first step to help guide, plan, 
and interpret a broader search for remnant radio galaxies with 
LOFAR. 

The paper is organized as follows: in Sect. [^ we describe the 
data and the data reduction procedures; in SectT|3]we present the 
source properties resulting from the data analysis; in Sect. we 
investigate the history of the source, addressing the question of 
how this object fits in the current scenario of active and remnant 
radio galaxies, and what its evolution history is. Finally, in Sect. 
I^we discuss the implications of this study for selecting remnant 
radio galaxies in upcoming surveys. The cosmology adopted 
throughout the paper assumes a flat universe and the following 
parameters: Hq = 70 km s 'Mpc ', Qa = 0.7, and = 0.3. 
The spectral index, a, is defined using the convention Fy - v““. 


2. Data 

2.1. LOFAR observations 


blobl has been detected with the High Band Antenna (HBA) 
within the field of view centred on the LOFAR primary calibra¬ 
tor 3C380 at a distance of about 30' from the phase centre. The 
observations took place on March 14, 2014 as part of cycle 1 for 
a total integration time of 8 h. The sampling time was set to 1 
second, and four polarization products (XX, YY, XY, and YX) 
were recorded. The total bandwidth, equal to 44.3 MHz, was 
divided by default into 227 sub-bands of 195.3 kHz with 64 fre¬ 
quency channels. Only the Dutch LOFAR stations were included 
in the observation, giving a maximum baseline of 100 km. We 
summarize the details of the observational configuration in Table 

[I] 


The data were pre-processed by the observatory pipeline 
( [Heald et’ak 2010 1 as follows. Each sub-band was automatically 
flagged for RFI using the AOFlagger ( [Offringa et al. 12012) 1 and 
was averaged in time to 5 seconds per sample and in frequency 
by a factor of 16, giving four channels per sub-band in the out¬ 
put data. The calibrator data were used to derive antenna gains, 
solving for a diagonal Jones matrix using the black board self¬ 
calibration tool (BBS, Pandey et al.|2009 1 , which takes into ac¬ 
count the variation of the LOFAR station beams with time and 
frequency. The source 3C380 was used as flux density calibrator, 
and the flux scale was set following Scaife & Heald ( 2012| l. 

The amplitude-corrected visibilities were merged in chunks 
of 20 sub-bands (3.9 MHz) and the channels were averaged to 1 
per sub-band. Each frequency chunk was imaged separately us¬ 
ing the awimager software (Tasse et al. 2013[ l and fixing the field 
of view to 5° X 5°. This imager makes use of the aw-projection 
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algorithm to model the primary beam variations with time and 
frequency across the field of view. Moreover, it corrects for the 
effects of non-coplanar baselines when imaging a large field of 
view. 

Phase self-calibration was performed iteratively, adding pro¬ 
gressively longer baselines to achieve the final resolution in three 
steps (Vilchez et al. in prep.). As a starting point, only core 
stations were used so that the phase coherence was preserved 
and because they are set to the same clock. The initial phase- 
calibration model was derived from the VLSS catalogue and 
contains all sources located in the field of view with flux den¬ 
sity higher than 4 Jy. For each VLSS source, counterparts in the 
WENSS and NVSS catalogues were searched for and associ¬ 
ated with these sources to provide spectral index information 
(see Sect. 6.4 in the LOFAR cookbool|^ for more details). We 
here took into account a first-order effect of the ionosphere in 
the direction of the dominant source, but did not perform any 
direction-dependent calibration. 

The images were produced using a robust weighting of -0.3 
and have a final beam size of 59 x 37". The final image was 
obtained by combining all images in the image plane, which 
resulted in a central frequency of 137 MHz and a nois e of 3.5 


mJy beam L For the spectral analysis presented in Sect. 3.4 


used two different sets of 20 sub-bands centred at 116 and 155 
MHz, respectively. The integrated flux densities of the source 
extracted from the maps are listed in Table The errors on the 
measurements were computed taking into account different con¬ 
tributions: i) 2% due to the uncertainty in the flux scale ( [Scaife 
|& Heald 2012|l, ii) 10% due to t he uncertainty in the global beam 
model ( |van Weeren et al.|2014| l, iii) 15% due to variations in the 
measured flux density when using apertures of different sizes, 
iv) ~ 10% due to the statistical error (thermal noise) measured 
around the source following [Klein et al.| ( |2003] l. 

Table 1. Details of the LOFAR observations. 


Obs. ID 
Date 

Phase centre 3c380 (J2000.0) 
Target coordinates (J2000.0) 
Total no. of stations used 
No. of remote stations 
Integration time 
Sampling time 
Polarization 
Frequency range 
Band width 
Synthesized beam 
(robust = -0.3) 

RMS (robust = -0.3) 


L206016 
14-March-2014 
18h 29m 31.8s, -H48d 44m 46s 
18h 28m 20.5s, -H49d 14m 20s 
60 
14 

8 hours 
1 second 
Full Stokes 
115-162 MHz 
44.3 MHz 

59 X 37" PA = -45° 

3.5 mJy beam^' 


2.2. WSRT observations 

Because blobl is not detected in the NRAO VLA Sky Survey 
(NVSS, Condon et al.|T998l l, we carried out deeper L4-GHz 
observations using the Westerbork Synthesis Radio Telescope 
(WSRT) on 12 and 24 August 2014. We also used the WSRT 
at high resolution to observe the field at 4.9 GHz on 29 March 
2015 to determine whether there are any compact components. 
The observations and data reduction are described below, and the 
setup is also summarized in Table 


Observations at 1.4 GHz 

In the first run, we observed for 12 h using the standard contin¬ 
uum setup of eight sub-bands of 20 MHz with 64 channels in 
each sub-band. Of these, four had technical problems and were 
discarded. In the second run, we observed for 10 h with six sub¬ 
bands of 20 MHz with 64 channels for each sub-band. An ex¬ 
tra sub-band of 20 MHz and 256 channels was dedicated to ex¬ 
plore possible Hi spectral lines and was centred on 1350 MHz 
to include HI at the redshift of the two candidate host galaxies 
(see Sect. 3.2). Primary calibrators (3C286 and 3C147) were ob¬ 
served at the beginning and end of each run. 

The data were calibrated and reduced using the MIRIAE|^ 
software package ( [Sault et al.|1995| l. Using standard cross cali¬ 
bration and self-calibration, the image is dominated by off-axis 
phase errors associated with a strong off-axis source (3C380), 
situated at ~ 30' from the field centre. These off-axis errors are 
most likely due to pointing or tracking errors and/or frequency- 
dependent errors in the primary model that are due to stand¬ 
ing waves in the dishes of the WSRT and cannot be corrected 
with standard direction-independent calibration. Therefore, we 
used the techni que commonly k nown as peeling, following the 


standard step^ (Noordam]2004 


to correct for off-axis errors in 
the image. The procedure involves the subtraction of the central 
sources in the field and self-calibration on individual problem¬ 
atic off-axis sources (in this case, one source, 3C380), iteratively 
producing antenna-based corrections and a source model for dif¬ 
ferent directions. The source model for 3C380 obtained in this 
way, including the additional off-axis corrections, was then sub¬ 
tracted from the original visibility data set. The peeling proce¬ 
dure was run separately on each 20 MHz sub-band. 

The final image of blobl was obtained in two ways. In the 
first case, we separately imaged each sub-band (i.e. 20 MHz) and 
combined the final images in the image plane. With the second 
approach, we produced an image by combining all the data in 
the uv plane. The two different strategies resulted in very similar 
maps. In both cases, we used a combination of robust weight¬ 
ing 0.4 and a 25" taper to achieve a low side-lobe level and a 
final resolution that allows imaging the low surface brightness 
emission expected for blobl. The final beam size is 30" x 27" 
(PA= -3°) and the noise of the image is 80 qJy beam '. The 
image is presented in Figj^ To image and locate compact com¬ 
ponents, we also produced an image with uniform weighting, 
reaching a resolution of 15" x 9" (PA= -4°) and noise equal to 
50 qJy beam '. The integrated flux density used for the spectral 
analysis was derived using only data from the first observing run, 
which have better quality. This is reported in Table The error 
on the measurements was computed taking into account differ¬ 
ent contributions: i) 2% due to the uncertainty in the flux scale, 
ii) 15% due to variations in the measured flux density when us¬ 
ing apertures of different sizes, iii) <1% due to the statistical er¬ 
ror (thermal noise) measured around the source following |Klein| 
|et al.| ( |2003| l, and iv) 20% due to variations in the measured flux 
density when using distinct cleaning strategies. 

No HI emission was detected. However, because subtracting 
3C380 was difficult, the H i cube has a high noise (about 0.8 mJy 
beam ' per channel with a channel width equal to 16.5 km s '), 
making the search for H i not sensitive enough for useful con¬ 
straints. The corresponding upper limit of Hi mass is derived 
following [Emonts et al.| ( [2010| l. At a redshift of z = 0.05 (see 
Sect|3.2|l and assuming a potential 3cr detection smoothed across 


' https://www.astron.nl/radio-observatory/lofar/ ^ 

http: //WWW .atnf.csiro.au/computing/software/miriad/ 

lofar-imaging-cookbook ^ 

http: //WWW .astron.nl/~oosterlo/peeling.pdf 
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a velocity range of 200 km s ' (typical for galaxies detected in 
H I), we derived an upper limit of Mhi < 1.5 X 10® Mq. 


Observations at 4.9 GHz 

The field was observed for 11 h in a standard continuum con¬ 
figuration (see Table 0. The array configuration was not opti¬ 
mized for detecting large-scale structure emission but for iden¬ 
tifying compact components at high resolution. A standard data 
reduction was performed using MIRIAD. The final image, ob¬ 
tained with uniform weighting, has a synthesized beam of 4.03 x 
2.94"(PA 0.5 deg) and a noise level of 85 juJy. 


Table 2. List of flux densities and respective errors. 


Telescope 

Frequency 

(MHz) 

Flux density 
(mJy) 

Error 

(mJy) 

LOFAR 

116 

1400 

300 

LOFAR 

155 

1200 

250 

WENSS 

327 

800 

150 

WSRT 

1400 

215 

45 

GB6 

4850 

<30 

- 


Table 3. Details of the WSRT observations. 


RA, Dec (J2000.0) 

18h 27m 0L9s -l-49d 12h 25.9s 

1.4 GHz - continuum 

Date 

12, 24-Aug-2014 

Integration time 

22 hours 

Frequency 

1400 MHz 

Band width 

4 X 20 MHz’-t- 6 X 20 MHz” 

Synthesized beam 


(robust=0.4) 

35 X 26" PA = -10.4° 

RMS (robust=0.4) 

0.08 mJy bearn^* 

Synthesized beam 


(uniform weighting) 

15 X 9" PA = -4° 

RMS (uniform weighting) 

0.05 mJy bearn^* 

1.4 GHz-H i 

Central frequency 

1350 MHz 

Band width 

20 MHz 

Number of channels 

256 

Velocity resolution 

16.5 km s^' 

4.9 GHz 

Date 

29 March 2015 

Integration time 

11 hours 

Frequency 

4900 MHz 

Band width 

8 X 20 MHz 

Synthesized beam 

4 X 3" PA = 0.5° 

(uniform weighting) 


RMS (uniform weighting) 

0.085 mJy beam^' 


Notes, ’usable bands in the first observation; ” usable bands in the 
second observation. 


2.3. Archival radio data 


We searched images of the relevant area from several radio sur¬ 
veys (i.e. VLA Low-Frequency Sky Survey, VLSS, [Cohen et aL 
( 2007| l; Faint Images of the Radio Sky at Twenty-cm Becker 
jet al. ( |1995| l and NVSS [Condon et al.j ( 1998[ l), and we only 
identified the source in the Westerbork Northern Sky Survey 
(WENSS, [Raigelink et al.|l997] ) at 327 MHz (see Fig.[2i. How¬ 
ever, we also obtained a useful upper limit at 4.8 GHz from the 
Green Bank 6 cm radio survey (GB6, [^ndon et al.|1994| l map, 
which was used in modelling the source spectrum (see Sect. 3.4 1 . 
The upper limit was extracted from the GB6 map and fixed at 3cr, 
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Table 4. List of two-point spectral indices and respective errors as a 
function of frequency. 


Frequency (MHz) 

Spectral index 

116-155 

0.50 ± 1.0 

155-327 

0.55 ± 0.40 

327-1400 

0.90 ± 0.20 

1400-4850 

> 1.60 


where cr was determined from a set of flux density measurements 
surrounding the source location using a box with the same size 
as the source. The parameters of these maps are listed in Table 

El 


3. Results 

3.1. Morphology 

The radio images obtained from the data described in Sect.|^are 
presented in Fig. [^ At all available frequencies, blobl shows a 
very diffuse amorphous shape emission with low surface bright¬ 
ness with a mean value of 4 mJy arcmin^^ at 1.4 GHz. Its size 
is 12' X 7', corresponding to an axial ratio of 1.7. We do not 
find evidence for edge brightening or sharp features such as jets 
or hot spots, but there is a broad linear enhancement of surface 
brightness extending to the north in both the 1.4-GHz WSRT and 
the LOFAR map. In the high-resolution 1.4-GHz WSRT image, 
which has a resolution of 15", we also identify a compact com¬ 
ponent with a flux density Si. 4 ghz = 1 mJy located near the cen¬ 
tre of the diffuse emission. This feature is also confirmed by the 
4.9-GHz WSRT observation with S 4 9 GHZ = 0.7 mJy (see Fig.[^. 


3.2. Optical identification and galactic environment 


Using data from the Digital Sky Survey (DSS), we identify 
two galaxies located near the centre of the diffuse radio source; 
these are the brightest optical objects within the field (see Fig. 
[^. The brightest of this pair of galaxies is co-located with the 
compact component detected in the WSRT images described in 
Sect. |3.1 and is considered the most likely host. The galaxies 


have been observed as part of the Two Micron All-Sky Survey 
(2MASS) and are catalogued as 2MASXJ18282048-H4914428 
(1) and 2MASXJ18282089-H4914278 (2). They have mi= 12.9 
and m2=13.8 in the Ks - band (Vega) at 2.17 pm, respectively. 

We obtained optical spectra for the two galaxies on 25 Febru¬ 
ary 2015 and deep imaging in a filter equivalent to the r band on 
30 April 2015 with the ISIS spectrograph and the ACAM instru¬ 
ment, both mounted on the 4.2 m William Herschel Telescope 
(WHT). The spectra were reduced using standard IRAlj^ tasks. 
The redshifts of the galaxies are zi - 0.051 and Z 2 = 0.052, 
both with relative uncertainties of Az = 0.00025, and absolute 
uncertainties of Az = 0.001. These are consistent with the pho¬ 
tometric redshifts from the 2MASS Photometric Redshift cata¬ 
logue (2MPR[Bilicki et al.|2014|l of Zphot.i = 0.053 + 0.006 and 
Zphot,2 = 0.050 + 0.006. 


The spectra of the two galaxies appear to be typical of early- 
type galaxies with strong 4000A breaks and stellar absorption 
features. No strong emission lines are detected in the spectra that 
might be connected to an optical AGN or to star formation. An 


IRAF is distributed by the National Optical Astronomical Observato¬ 
ries, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Sci¬ 
ence Foundation 
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29m00.00s 40.00s 20.00s 28m00.00s 18h27m40.00s 

RA 02000) 



Fig. 1. Radio maps of blobl: (top panel) LOFAR 137-MHz map, levels 
2, 3, 5, 10 X o- (3.5 mJy beam-*); (middle panel) WENSS 327-MHz 
map, levels 2, 3, 5, 10 x cr (5 mJy beam-*); (bottom panel) WSRT 
1400-MHz map, levels: 2, 3, 5, 10 x cr (0.08 mJy beam-*). 


derived from this portion of the spectrum. The HjS- [OIII] region 
in either spectrum shows no signs of AGN activity. 

If the radio emission is associated with the two galaxies, it 
would imply a linear size for the radio emission of 700x410 kpc. 
In Fig. I^the 1.4-GHz WSRT radio contours are overlaid on the 
2MASS map, showing the extent of the radio emission. From the 
Ks - band luminosity and assuming a mass-to-light ratio of 0.6 
( |McGaugh & Schombert|[2014| l, we derived a hrst-order stellar 
mass estimate of logMi = 10.6 Mq and logM 2 = 10.3 Mq. 

The close physical separation and radial velocity difference 
of ~ 300 km s“' suggest that the galaxies are interacting. This is 
conhrmed by tidal features observed in the r-band image, as is 
shown in Fig.j^and further discussed in Sect. |4.2| 

The two galaxies do not belong to any known galaxy cluster 
or group. The closest galaxy cluster is MCXC J1811.0+4954 
(z=0.0501), which is located at a distance of 2.7° 10 Mpc 

( |B6hringer et al.|2000) l. In our deep r-band image we identify 
at most three or four candidate members within a projected ra¬ 
dius of 50 kpc from the host galaxy. No spectral information is 
available for these objects, therefore we cannot verify if they are 
background or foreground sources. We therefore conclude that 
the pair of galaxies might be part of a small group, but probably 
not part of a massive cluster. 

3.3. Source energetics 

In this and the following section, we derive physical parameters 
of the radio emission and model the radio spectrum to assess 
the nature of blobl and how it compares with other sources. 
Assuming J18282048+4914428 to be the host galaxy, the ra¬ 
dio luminosity of the diffuse emission at 1.4 GHz is Li 4 ghz = 
1.5 X 10^4 WHz-*. 

We computed the magnetic held strength to be 1 gG and the 
total energy density to be 1.6 x 10“'^ erg cm“^ , assuming the 
plasma to be in equipartition condition between particles and the 
magnetic held using the derivation of Beck & Krause ( 2005) 1. 
This implies a non-thermal pressure of 5 x lO"^'* dyne cm For 
this calculation we assumed a power-law particle distribution of 
the form N(y) cx y-** between a minimum and maximum Lorentz 
factor of ymin = 10 and y^ax = 10®. We calculated the volume 
of the source as 1.7 x 10^^ cm^ assuming a prolate ellipsoidal 
geometry with the major and minor axis equal to the major and 
minor projected size of the source. We furthermore assumed the 
particle energy content of the source to be equally distributed 
between heavy particles and electrons so that their ratio k=l, 
and we set a = 0.5 according to the observed low-frequency 
spectral shape where radiative energy losses are less signihcant. 
We note that with this low value for the magnetic held, the ra¬ 
diative cooling of the plasma is dominated by inverse Compton 
(IC) scattering of cosmic microwave background (CMB) pho¬ 
tons: the energy density of the CMB is a factor 10 greater than 
the magnetic energy density. 


upper limit on the star formation rate can be derived from the 
[Oil] emission line Hux, using the reddening-corrected relation 
between [Oil] luminosity, L[OII], and star formation rate derived 
in Eq. 4 from Kewley et al. ( 2004jl . A synthetic unresolved [Oil] 
3727,3730A doublet of width 3A and peak Hux equal to 3-cr 
of the noise measured between 3800-4100A (3615-3900A rest- 
frame) was inserted into the spectra of the two galaxies. The star 
formation rates resulting from this emission line are SFRi = 0.51 
and SFR 2 = 0.15 Mo/year. Unfortunately, Ha falls directly on a 
telluric absorption line, which means that no information can be 


3.4. Spectral properties and age 

The spectral shape of radio sources is directly related to the elec¬ 
tron populations responsible for the emission and gives impor¬ 
tant constraints about their age, as we describe below ([Keller- 1 


1964) Pacholczyk 


1970 1 . The integrated radio spectrum 


of active galaxies is well described by the continuous injection 
model (Cl, |Kardash^| 1 962} , which assumes a continuous re¬ 
plenishment of new relativistic particles with an energy distribu¬ 
tion such that N(E) oc E **. This results in a broken power-law 
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radiation spectrum with spectral index ainj = ^ below a criti¬ 
cal frequency Vb,iow and a - ainj + 0.5 above Vb.iow that is due to 
preferential radiative cooling of high-energy particles. The typi¬ 
cal ainj of active sources is in the range of 0.5-0 .8 and, as such, 
has spectral index a = 1 ~ 1.3 at high frequency. When the nu¬ 
clear activity switches off, a new break frequenc y Vb.high appears 
beyond which the spectrum drops exponentially ( [Komissarov &| 
Gubanov|1994| l. Knowledge of the spectral shape and the source 
energetics allows estimating the age of the electron population 
according to radiative cooling as follows: 

^ 0.5 

h = 1590—- " V , (1) 

(Bl^ + Bcmb) Vvb(l +z) 

where f^ is given in Myr, the magnetic field B and inverse Comp¬ 
ton equivalent field Bcmb are in pG, and the break frequency Vb 
is given in GHz. Furthermore, the ratio between Vb.iow and Vb^high 
constrains the active time of the source with respect to the total 
age as follows: 


tOFF 


t^b.high 

^b.low 


0.5 


( 2 ) 


where tj is the total age of the source and can be written as tj = 
tci + toFF , where tci is the continuous injection phase and toFF is 
the off phase. 


As shown in Fig. ^ the spectrum of blob 1 is highly curved 
at a frequency > 1.4 GHz. The measured spectral indices at 
low frequency are in accordance with typical observations of 
active radio sources (see Table |^. At higher frequencies, the 
spectral shape steepens to reach a value of aj^QQ > 1.6. Fol- 
lowing jMurgia et ^pOl Ij l, we computed the spectral curvature 
SPC = ahigh - aiow using aiow = a^l and aygh = a|^^° and 
obtained a value of SPC > 1.1. This spectral curvature is well 
beyond the values expected for active radio sources according 
to the continuous injection model discussed above (SPC < 0.5), 
therefore it suggests that the source of the particle injection is 
switched off. 

Because the slope variation of the spectrum is so abrupt, we 
can assert that the two break frequencies Vb.iow and Vb.high He in 
this frequency range and are close to each other. This implies, 
according to Eq. that the source has been inactive for a sig¬ 
nificant fraction of its total age. Assuming viow >1-4 GHz and 
Vhigh ^ 5 GHz, we can calculate to first order the total age and 
the off phase of blob 1 using Eq[T]for synchrotron cooling. From 
this we derive the following limits for the source total age and 
off phase: ts < 95 Myr and Ioff ^ 45 Myr. 

To obtain a more accurate estimate of the activity timescales 
of the source, we fit the spectrum with a model describing the 
spectral evolution of a source that is active for a time tci and then 
switches off for a time toFF = ts “ tci ([Komissarov & Gubanov| 


1994| Murgia et al. 2011 1 . For the modelling we assumed a fixed 
injection index of ainj = 0.5 and the magnetic field value 1 pG 
derived in Sect. |3.3| The best fit is shown in Fig.j^and the best- 


fit parameters are toFF = 60 ;^ 2 o Myrs and ts/toFF = l- 3 ;[;g 3 (see 
Table|^. This implies an active phase of tci = a 

total age of ts = 75Myr, which agrees with the previous first- 
order estimates. We note that because of the assumptions and 
systematic errors that enter the equation, these numbers should 
only be considered as characteristic timescales. 
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Table 5. Best-fit parameters of the spectral modelling presented in Sect. 

nn 


ffmj B 

ts/toFF 

tOFF 

pG 


Myr 

0.5 fixed 1 

1 

^^-20 



Frequency (MHz) 

Fig. 3. Integrated radio spectrum of blobl. The solid line represents 
the best fit of the synchrotron model described in Sect. |3.4| The flux 
densities used for the fitting are summarized in Tablej^ 


4. Discussion 


4.1. Evidence in favour of the remnant radio gaiaxy 
interpretation 

The morphology of blobl with its amorphous shape and low- 
surface brightness together with its spectral curvature supports 
an interpretation in which the source is a remnant of a past radio 
AGN activity. Here we discuss these aspects in detail in the con¬ 
text of current knowledge of the population of active and rem¬ 
nant radio galaxies. 

According to the P-D (power-size) diagram for low-power 
radio galaxies by jParma et al.| ( [T999| l, blobl has a relatively large 
size for its luminosity. This may suggest that the central engine 
was more powerful during the maximum of its activity. As the 
jets switch off, powerful radio galaxies are expected to undergo 
a rapid luminosity decrease (jSlee et al.||2001| Kaiser & Cotter 


2002| Hardcastle & Krause 2013| l. Another explanation could 
be that the radio source has expanded in an environment with a 
comparatively low density. However, blobl still lies within the 
broad scatter of the P-D correlation, therefore this is not conclu¬ 
sive. 

The core of blobl observed at 1.4 GHz and 4.9 GHz has a 
power of 5.7x10^' W Hz ' at 1.4 GHz, giving a core prominence 
(i.e. ratio between the radio power of the core at 5 GHz and the 
extended emission at 178 MHz) of about 6 x 10 This value 
lies more than a factor 10 below the correlation observed for 
active radio galaxies by [Giovannini et ar] ( |1988| ). This suggests 
that the nuclear engine may be still active in a low-power state 
and that the broad linear feature presented in Sect. |3.1| might be 
related to an uncollimated weak outflow connected to the core. 
The spectral index of the core is = 0.3 in accordance 

with its AGN nature. Star formation is unlikely since the 1.4 GHz 
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core lum inosity would transla te into a star formation rate of 1.2 
Mq yr ' ( [Condon et al. 2002 1 , more than twice as much as what 
is derived from the optical spectrum. The lack of strong AGN 
emission lines is also consistent with the central engine being 


in a dying phase or in a low-power ’hot mode’ state (Tadhunter 
let al.|1998» . 

Remnant radio galaxies are expected to lack compact fea¬ 
tures typical of active galaxies, such as defined jets or hot spots. 
Depending on whether the central activity has ceased gradually 
or abruptly, they can also show weak or absent cores. The mor¬ 
phology and evolution of the remnant does depend on the exter¬ 
nal medium and on the properties of the progenitor active source. 


which are different for FRI and FRII radio galaxies. Saripalli 
|et al.| (|2012 1 suggested the following criteria for distinguishing 


et al. 


( 2007|l showed that such fat bubbles of plasma can be in¬ 


flated by wide and slow jets. It is worthwhile mentioning that 
ISubrahmanyan et ah] ( |2006| ) also claimed that the relaxed or ex¬ 
tended morphology observed in some sources might be the re¬ 
sult of jet precession during the fly-by encounters that precede a 
merger event. In this case, the synchrotron plasma is deposited 
over a wide range of angles conferring the observed morphology. 
As discussed in Sect. 3.2 the host galaxy of blobl is interacting 
with a companion, therefore this last scenario may also be rele¬ 
vant to the evolutionary history of blobl. 

4.2. Merger scenario and triggering mechanism 

A complete investigation of the evolutionary history of blobl 
involves the analysis of the mechanism responsible for trigger¬ 
ing and stopping the radio-loud AGN. Because the host galaxy 
of blobl is currently interacting with another galaxy, we should 
consider the possibility that this process had a role in marking 
the timescales of the radio activity of the AGN. 

Galaxy mergers are often suggested to influence the trigger¬ 
ing and quenching of the jets in powerful radio-loud AGN ( jTad-j 
hunter et al.|2012| Wen et al.|2012 1 . They allow the gas to lose 
angular momentum and flow towards the nuclear region of the 
galaxy, whereby they either contribute to replenishing the black 
hole fuel, which in turn enhances the jet production, or to unset¬ 
tling the accretion, which causes the disruption of the jets. Tidal 
encounters between galaxies that precede the final coalescence 
can also be responsible for triggering or stopping the AGN ac¬ 
tivity ( |Ellison et al.|2011~| l. In this phase, intermittent jets can be 
produced, and if the accretion disk is perturbed, they can also be 
realigned. One of the best-studied example of this occurrence is 
CenA, where a merger is thought to be responsible for the inter¬ 
ruption of the radio activity for a time of the order of 10 Myr, af¬ 
ter which new-born jets have restarted ( Saxton et al.|200lj Struve 
|et al.||20T0|l. Although statistical evidence of the link between 


FRI from FRII remnant sources on a morphological basis: i) 
FRII remnant sources appear as pairs of diffuse edge-brightened 
lobes lacking hot spots, jets, and cores; ii) FRI remnant sources 
do not show any compact emission except for a weak core with 
twin extensions or trails that do not resemble typical jets or lobes. 
With its relaxed lobes (axial ratio of 1.7), absence of edge bright¬ 
ening, and weak core, blobl is closer to an FRI morphology. 
According to Parma et al. ( |1996| l, about 62% of the B2 sam¬ 
ple of low-luminosity radio galaxies have double-lobe morphol¬ 
ogy with relaxed lobes whose typical axial ratio is 2 , whereas 
high-power FRII radio galaxies show much higher axial ratios 
of about 5. The standard prototypes of the “fat double” sources 
(Owen & Laing 1989) are identified with 3C310; 3C314.1 and 
3C386 having power in the ranges 178 = 10^"^-10^^ W Hz 'sr ' 
and axial ratios < 2. Hydrodynamical simulations by jSternber^ 


mergers and radio emission has been found, the AGN radio ac¬ 
tivity has not been observed to be associated to a specific phase 
of the interaction ( |Ramos Almeida et al.|201 1) . 

From the r-band image it is clear that the morphology of 
the two galaxies is disturbed. The host galaxy shows an asym¬ 
metry in the brightness distribution on the top east side while 
the companion shows clear tidal features on the west and east 
side. Although mergers with late-type galaxies can produce the 
most prominent tidal features, early-type galaxies are also ob¬ 
served to produce tidal tails ( Ramos Almeida et al.||20lT] l, but 
much broader, diffuse, and shorter lived than what is observed in 
gas-rich mergers. Since both galaxies are early-type, we do not 
expect much gas to be involved in the interaction and therefore 
not much star formation to be going on. This agrees with the 
computed star formation rate presented in Sect. |3.2| According 
to their morphology, it is likely that the two galaxies are observed 
after the first pericentre passage but before the final coalescence 
of the merging nuclei. 

We estimate the dynamical timescale of the system t^yn = 
60 Myr based on the ratio between the projected separation be¬ 
tween the two galaxies dpioj = 15 kpc and their relative speed 
Vrei = 300 km s'. Although this timescale is uncertain and needs 
to be taken with care, its similarity to the topF of the radio source 
may suggest that the interaction could have disrupted the radio 
jets. 


4.3. Spectrai comparison with other remnant radio gaiaxies 
and impiications for remnant evoiution 


Based on the discussion in Sect. |4.1| we consider blobl to be 
a remnant radio galaxy and compare here its spectral character¬ 
istics with those observed in previously selected sources of the 
same class. In particular, we mainly focus our comparison on 
the only complete samples of dying radio galaxies reported by 
jParma et al. ( |2007| l and jMurgia et al.j ( |201 Ij), for which a uni¬ 
form multi-frequency analysis is available. Parma et al.|(|20071l 
selected their sample using 0 : 327 *^ > 1.3, while Murgia et al. 
( | 2011 ) l selected their sample based on the spectral curvature de- 
flned in Sectl3.4l 


Although it is generally believed that remnant radio galaxies 
develop spectral indices > 1 below 1.4 GHz, this is not the rule, 
and blobl is one of the few known examples of remnant sources 
with low-frequency spectra consistent with active sources. A 
similar case is the source B2 1610-1-29 in Murgia et al. ( 2011| l. 
These remnant galaxies show the high-frequency spectral cut-of 
that is typical of the dying phase only above 1 GHz, while at 
lower frequency they still preserve the injection spectral index 
that characterizes the active phase. This spectral shape can ei¬ 
ther be justified by the plasma being still young or by a slow rate 
of radiative cooling. 

The spectral modelling of blobl presented in Sect. |3.4| shows 
that the source has been spending about 50-80% of its life in a 
dying phase. The computed total age agrees with the mean value 
observed by |Murgia et al.| ( |201 l| l and jParma et al.| ( |2007[ ) and 
also with what has been derived in other individual sources such 
as B2 0924-1-30 (60 Myr; jJamrozy et al.||2004| Shulevski et al. 
in prep), B 0917-H75 (100 M yr; |Harris et al.|1993l l, J1324-3138 
(97 Myr; Venturi et al.|1998 1 , and A 2622 (150 Myr; Giacintucci 
|et al.|2007| l. The same can be asserted for the ratio topp/ts, which 
fits the range of values 0.2-0 .8 computed in Murgia et al. ( 201 l) l 
and Parma et al. ( 2007| l within the errors. The difference in spec¬ 
tral shape can thus be explained by different values for magnetic 
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fields (sources in Murgia et al. (2011 1 have Beq =4-45 pG 
compared with 1 fiG in blobl). 

It is interesting to note that despite being associated with a 
field galaxy, blobl has a long dying phase whose duration is 
comparable to those observed for remnants in clusters. Dense 
environments are thought to keep the plasma confined, prevent¬ 
ing adiabatic expansion and allowing their detection for a longer 
time. Instead, in low-density environments, the plasma is free to 
expand, and if it is over-pressured with respect to the external 
medium, the source can rapidly fade away. This argument is of¬ 
ten invoked to explain why so few remnants are detected outside 
clusters, although a difference in the duty cycle of radio galaxies 
in different environments can also be an explanation (e.g. |Mur^ 
|gia et al.]|2011| l. To date, very few other remnant sources out¬ 
side clusters have been detected and studied; 0924-1-30 Pamrozy 


let al.|2004|), NGC 5580 (|de Gasperin et al.|2014| l, and NGC 1534 
(|Hurley-Walker et al.|2()15|[ 


The reason why the plasma is still visible after such a long 
time has to be connected with its energetics. If the lobes have low 
internal energy densities and pressures, their expansion rate is 
low and the radiative cooling is the dominant process in the evo¬ 
lution of the source. This is very likely the case for sources with 
a low surface brightness like blobl ( {Subrahmanyan et al.|2003| l. 
Furthermore, it has been argued that most sources will be at pres¬ 
sure equilibrium at the end of their lives ( |Hardcastle & Krausej 
|2013^ , but without further information on the external medium 
we cannot evaluate whether the plasma is over-pressured or in 
pressure balance. The only archival X-ray observations of this 
area of sky are from the ROSAT survey, and because they are 
too shallow, they do not yield useful constraints. 

Better comprehension of the phenomena that drive the 
plasma evolution as a function of energy and environment can 
only be achieved with larger samples of this class of sources. For 
this reason new systematic searches exploiting new-generation 
instruments such as LOFAR need to be performed, as we dis¬ 
cuss in more detail in the following section. 


5. Implications for selection of remnant radio 
galaxies 

Until now, systematic searches of remnant radio galaxies have 
mostly been carried out by cross-matching radio surveys at dif¬ 
ferent frequencies using the steep spectral index as source selec¬ 
tion criterion (typically a > 1.3). |Parma et al.| (2007 1 selected 
a sample of dying radio galaxies from 327-MHz WENSS and 
the 1.4-GHz NVSS catalogues. Dwarakanath & Kale ( 2009) 1 and 
van Weeren et al.|(|2009)l extracted a sample of remnant candi¬ 
dates from the 74-MHz VLSS and 1400-MHz NVSS surveys. 
A few tentative searches have also been performed on individ¬ 
ual deep fields (Sirothia et al. 2009 [ |van Weeren et al. 2014| l. 
The strategy of the steep spectral index is best for identifying 
sources with very curved spectra at low frequencies that have 
experienced significant radiative losses and are typically not de¬ 
tectable at 1.4 GHz. 

The remnant radio galaxy presented in this paper demon¬ 
strates that depending on its physical conditions and environ¬ 
ment, the plasma can undergo a different evolutionary path and 
present different spectral shapes. Therefore a number of sources 
might be missed by this standard selection approach and new 
selection methods should be considered. |Murgia et al.| ( |2011| l 
suggested that the SPC can be a powerful diagnostic tool for 
tracing the evolutionary stage of radio sources. On the basis of 
synchrotron models, dying sources are expected to have an SPC 


> 0.5. With this method, they selected a sample of five dying ra¬ 
dio galaxies by combining multi-frequency survey data and dedi¬ 
cated observations. The SPC measured for blobl supports the re¬ 
liability of the strategy. The broad range of frequencies required 
to compute it makes it hard to apply it systematically as selection 
criterion, however. The lack of deep, all-sky, high-frequency sur¬ 
veys (~5 GHz) represent a limitation for this kind of approach. A 
combination of low frequency (e.g. 60-600 MHz) and high fre¬ 
quency (e.g. 1400-8000 MHz) is indeed needed to properly trace 
the spectral curvature ( |Harwood et al.|[2()13| l. We also note that 
sources with power-law non-steep spectra like the one studied 
by de Gasperin et al. ( 2014| l would likewise be missed with this 
method. 

In light of this, we suggest that the population of remnant 
radio galaxies cannot be fully characterized by a single selec¬ 
tion criterion based on spectral index or spectral curvature, for 
instance, but these should rather be used in a complementary 
way. Given the variety of observed spectral shapes and to con¬ 
struct complete samples, morphological selection should also 
play a role alongside spectral criteria. Visual inspection of radio 
sources has already provided a number of detections in addition 


to the source presented here (e.g. 

de Gasperin et al. 2014 

{Harris 

|et al.|1993| {Hurley- Walker et al. 

2015[ Giovannini et al.|1991|l. 


but it still has not been apphed in a systematic way. Tentative 
searches based on morphological characteristics have only been 
performed by |Saripalli et al.| (|2012|) in the context of the Aus 


tralia Telescope Low Brightness Survey and by Jones & Pre 


ston 


(20011 using the VLSS. Because remnant sources are ex¬ 


pected to be very diffuse with very low surface brightness emis¬ 
sion, new-generation instruments such as LOLAR with unprece¬ 
dented sensitivity capabilities are expected to unveil a new pop¬ 
ulation of such sources. The mean surface brightness of blobl 
(equal to 4 mJy arcmin^^ at 1.4 GHz) is indeed among the low¬ 
est values ever observed, and it lies just at the border of what can 
be detected with currently available surveys. This value is only 


33 (3 mJy arcmin 


Subrahmanyan et al. 2003 

1 , the remnant 

/ arcmin ^ 

Shulevski et al. 

2012 

1, several 


restarted LR-I and LR-II radio galaxies detected in the Aus tralia 
Telescope Low Brightness Survey (2.5-6 mJy arcmin^^, Sari- 


Ipalli et al.||2012| l, and a few ac tive giant radio galax ies such as 
J0034.0-6639 (1 mJy arcmin^^ 

J0515-8100 (10 mJy arcmin-^, 


, Sarinalli et al. 

2012 

1 and SGRS 

Subrahmanyan et al. 

2006 1. 


In conclusion, deep imaging at low frequency will surely 
probe a fainter population of steep-spectrum sources, which will 
expand the previous samples obtained from VLSS, WENSS, and 
NVSS, but will most likely also reveal diffuse sources with a flat¬ 
ter spectrum, like blob 1. To avoid selection biases towards steep 
integrated radio spectra, additional selection criteria should play 
a role in identifying remnant AGN candidates; these would be 
spectral curvature and morphology. A complete representation 
of the remnant radio galaxy population can be provided only in 
this way. By increasing the detection of these rare sources, we 
will finally be able to perform uniform systematic studies that 
will improve our knowledge of the different evolution histories 
of radio galaxies and of the different mechanisms that drive the 
remnant plasma evolution after the jets switch off. 
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Fig. 2. Left - blobl contours from the 1.4 GHz WSRT map (levels 3, 6, 12, and 15 tr) overlaid on the 2MASS Ks - band image of the field; right - 
Zoom-in on the host galaxy and on the companion galaxy from the r-band image contours from the 5 GHz WSRT map in black (levels 3, 5, and 7 
X <r (0.085 mjy beam^*). 
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